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A localized qubit entangled with a propagating quantum field is well suited to study non-local 
aspects of quantum mechanics and may also provide a channel to communicate between spatially 
separated nodes in a quantum network. Here, we report the on demand generation and character- 
ization of Bell-type entangled states between a superconducting qubit and propagating microwave 
fields composed of zero, one and two-photon Fock states. Using low noise linear amplification and 
efficient data acquisition we extract all relevant correlations between the qubit and the photon states 
and demonstrate entanglement with high fidelity. 



One of the most fascinating aspects of quantum physics 
is the entanglement between two spatially separated ob- 
jects sharing a common non-local wave function. Prop- 
agating photons are ideal carriers for distributing such 
entanglement between distant matter systems in a quan- 
tum network. Entanglement between photons and sta- 
tionary qubits has so far been exclusively studied at op- 
tical frequencies with single atoms [1-3 and electron 
spins [4 , to interface stationary and flying qubits [4, 
to implement quantum teleportation [6l |7] and to real- 
ize nodes for quantum repeaters [8 and networks . 
Rapid progress in the development of superconducting 
circuit based quantum technologies also renders propa- 
gating p!2UT6] and localized microwave photons [171 [IH] 
an attractive carrier of quantum information. A major 
obstacle in measuring quantum correlations between su- 
perconducting artificial atoms and itinerant photons has 
so far been the limited detection efficiency at microwave 
frequencies. Here, we overcome this problem by using 
a quantum limited amplifier [19j, which significantly im- 
proves the signal-to-noise ratio in both photon field and 
qubit measurement. In combination with novel tomog- 
raphy methods [20 this allows us to measure quantum 
correlations between itinerant microwave radiation and a 
stationary qubit with high fidelity. 

In our experiments we create entangled states between 
a superconducting transmon qubit [21] and an itinerant 
microwave radiation field containing up to two photons. 
The field mode a can be described by two canonically con- 
jugate variables X and P analogous to the position and 
momentum variables of a mechanical quantum harmonic 
oscillator. In contrast to most experiments performed at 
optical frequencies, we simultaneously measure both con- 
tinuous variables X, P rather than the photon number of 
the field. This enables us to fully characterize quantum 
fields beyond the single photon level. In addition to the 
measurement of photon statistics of the field, we deter- 
mine the correlations between the measured qubit state 
and the observed values X and P which clearly demon- 
strate that the qubit is entangled with the quantum field. 

We deterministically prepare 1.25 X 10^ Bell states of 
the form = (|0e) + \lg))/V2 per second, in which 
a single excitation is shared coherently between a qubit 
and a single propagating mode of a radiation field. Here, 
1^), |e) label the qubit basis states and |0), |1), |2), ... 
the photon number states. To entangle the qubit and 
the radiation field we first bring the qubit from the 



ground state |0^) to the excited state |0e) by applying a 
10 ns long TT-pulse resonant with its transition frequency 
uOge/27: = 6.442 GHz. By applying a magnetic flux pulse 
we then tune the qubit into resonance with a trans- 
mission line resonator at frequency co'r/27r = 7.133 GHz, 
which is strongly coupled to the qubit with rate g/27r = 
65 MHz [22]. After an interaction time of r = 7r/4^ ^ 2 ns 
we obtain the state (|0e) + \ lg))/\^ up to a phase factor 
which is omitted here for convenience. 

The setup employed for this experiment is shown 
schematically in Fig. [ija) with a micrograph of the sam- 
ple shown in B and the experimental sequence in D. 
While the qubit on average keeps its excitation during its 
life time Ti = 1.0 /is, the resonator field is emitted during 
the much shorter cavity decay time of 1/a>: = 25 ns. Note 
that, the interaction time r during state preparation is 
small compared to 1//^:, which itself is small compared to 
the qubit life and coherence times (T2* = 220 ns). This 
hierarchy of timescales {1/g < < Ti,T2*) guaran- 
tees that the entangled state can be coherently prepared, 
and that the qubit remains in the excited state while the 
photon is emitted into the propagating transmission line 
mode a of which both conjugate field quadratures X and 
P are detected (Fig.jljd)). 

We measure the quadratures X, P using a parametric 
amplifier operating close to the quantum limit [19 . The 
amplifier is based on a quarter wave transmission line 
resonator shunted by an array of Josephson junctions 
providing the Kerr nonlinearity used in the parametric 
amplification process |23], see Fig. [TJc). We operate the 
parametric amplifier in a phase-preserving mode in which 
both conjugate field quadratures are amplified equally 
[24] . This is achieved by pumping the amplifier at a fre- 
quency 12.5 MHz detuned from the center frequency of 
the detected photon pulse uj^l2i\. Note that we have op- 
erated the parametric amplifier with a relatively moder- 
ate gain of G(cJr) = 16.5 dB at the center frequency of the 
radiation field to be detected, in order to have a flat gain 
curve over the entire band of detection. In this regime, 
the effective detection efficiency is an order of magni- 
tude higher than for typical setups using transistor based 
amplifiers only After amplification, we record the 

time dependent quadrature amplitudes {X(t\P(ty\ in a 
microwave frequency heterodyne detection setup similar 
to the ones discussed in Ref. [22^. Performing temporal 
mode matching by convolving {X(t),P(t)} with an ap- 
propriate filter function [W, we retain one pair of values 
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{X, P} per generated Bell state. 

After the field detection we perform qubit state tomog- 
raphy (Fig. [lid)). We measure the qubit Bloch vector 
components ((Ja^), (ay) and (a^) by rotating the qubit 
into the respective eigenbasis and then applying a co- 
herent read-out tone to the resonator. Due to the dis- 
persive resonator frequency shift of x/^tt = 2.1MHz the 
integrated phase quadrature Q of the transmitted time- 
dependent signal Q{t) depends on the measured qubit 
state [25H27] . The probability distribution p{Q) is fitted 
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FIG. 1: Schematic of the experimental setup, (a) A trans- 
mon qubit with individual charge drive line (red) at frequency 
(jjge and flux control line (green) is strongly coupled to a res- 
onator with a weakly coupled input port (7in/^ ~ 0.001) 
driven at frequency cjr for qubit readout (orange) . The output 
port is strongly coupled into a transmission line (violet) and 
amplified with a Josephson parametric amplifier (light blue) 
pumped at ujp through a directional coupler with adjustable 
phase and attenuation for pump tone cancelation. The signal 
reflected off the parametric amplifier passes through a chain 
of circulators into a low-noise semiconductor-amplifier after 
which its two quadratures X, P are detected, (b) False color 
optical micrograph of the sample. The transmon qubit (en- 
larged) consists of two capacitively coupled islands connected 
by a pair of Josephson junctions, (c) False color micrograph 
of the Josephson parametric amplifier. The array of Joseph- 
son junctions (enlarged) at the end of the quarter wavelength 
resonator (light blue) provides the nonlinearity. (d) Pulse se- 
quence used for the experiment: (i) state preparation, (ii) field 
measurement and (iii) qubit read-out (see text for details), (e) 
Measured probability distribution p{Q) of the qubit read-out 
quadrature Q for prepared ground (blue), excited (red) and 
Bell (dashed white) states. 



to a weighted sum of two independently measured ref- 
erence distributions for the ground and excited state to 
extract the excited state population in the chosen basis 
(Fig. [TJe)). Due to qubit decay during the time required 
for the measurement of photon field quadratures after 
preparation of the entangled state, the single-shot qubit 
readout fidelity making use of the same mode is limited 
to 37%. In future experiments this aspect could be im- 
proved by using separate modes for photon generation 
and qubit read-out, similar to Ref. [28 . 

We extract the correlations between qubit and photon 
in the generated Bell states by recording 3-dimensional 
histograms of triplets {X, P, Q}, which count the num- 
ber of times for which the qubit read-out quadrature Q 
is measured in combination with photon field quadra- 
tures X and P. Within the limitations of the avail- 
able memory, we chose to discretize the histograms into 
128 X 128 bins for the measured photon field quadratures 
times 8 bins for the measured qubit quadrature. Efficient 
data acquisition and generation of histograms is realized 
in real-time using field programmable gate array elec- 
tronics [13 . The resulting data, which is obtained after 
preparation and detection of ~ 3 x 10^ Bell states, con- 
tains complete information about the photon statistics 
as well as all relevant qubit-photon correlations. From 
the measured histograms we extract the qubit population 
for each quadrature pair {X,P} by fitting the histogram 
columns along the Q axis to the two reference histograms 
shown in Fig. [TJe). 

Preparing two reference states |0^) and |0e) - for which 
the photon field is left in the vacuum state and thus is 
not correlated with the qubit - we find that the qubit 
population (blue: ground state, red: excited state) is 
independent of the detected field quadratures X and P 
(see Fig. [2|a) (i)-(ii)). This also indicates that there are 
no correlations of technical origin between the qubit and 
photon field measurements. The standard deviation of 
the photon field distribution 5m = 1.84 is larger than the 
quantum limit 1/a/2, due to noise added by the ampli- 
fiers, losses in the cables and microwave components, as 
well as finite mode matching efficiency, the combination 
of which corresponds to an effective detection efficiency 
of T] = 15%. Individual measurement results with large 
amplitude values (V-^^ + P^ ^ 3(5^) are unlikely, which 
causes the larger statistical uncertainty in the extracted 
qubit populations at the boundary of the colored regions. 

When Beh states |?/^) = (|0e) + \lg))/V2 are pre- 
pared we find a clear dependence of the measured qubit 
Bloch vector on the measured field quadratures {X, P} 
(Fig. [2|a) (iii)-(v)), in stark contrast to the results ob- 
tained for separable states. Measuring the qubit in the 
az basis we find a higher probability to observe the qubit 
in its ground state at large measured field amplitudes 
(blue region) and a higher probability to find the qubit 
in the excited state at small measured field amplitudes 
(red region in Fig. [2](a) (iii)). This observation is consis- 
tent with the expectation to either find the qubit in the 
ground state when a photon is propagating in mode a or 
in the excited state when a is in the vacuum state result- 
ing in correspondingly small field amplitudes. The fact 
that the measured qubit population is circularly symmet- 
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FIG. 2: Photon-qubit correlations for a prepared Bell state, (a) Qubit state population conditioned on the measured photon 
field quadratures X and P for the indicated Bloch vector components. {X, P} pairs for which no measurement results occured 
are shown in gray. White circles indicate the standard deviation of the photon field distribution 8^- (i)-(iii) Wz) for the 
reference states lO^^), |0e) and the Bell state IV^). For better visibility only, the data in subpanel (iii) is offset by its total 
mean ^ 0.1 compensating for the qubit decay during photon detection, (iv)-(v) {<Jx) and {ay) for the Bell state l^^). (b) 
Expectation values {{a^)^a^ai) extracted from the qubit-photon field correlations shown in A and the measured photon field 
distribution. The real (imaginary) part of these measured moments are shown in red (blue) and compared to the ideal Bell 
state (wireframe). The error bars are extracted from the standard deviation of repeated measurements, (c) Real part of the 
measured (solid) and ideal (wireframe) density matrix p for the Bell state l-^) with fidelity F — = 83%. 



ric in phase space, i.e. it is independent of the phase of 
the propagating field, indicates that \g) and |e) are cor- 
related with Fock states - such as the single photon or 
the vacuum state. 

To distinguish the coherent superposition of |0e) and 
\lg) in the Bell state from a mere statistical mixture, 
we measure the equatorial components {(Jx) and {(jy) of 
the qubit Bloch vector by applying 7r/2 pulses to the 
qubit about the corresponding axes and determine their 
correlations with the measured radiation field. We find 
that whenever a positive field quadrature X > is mea- 
sured the qubit is more likely to be found in a state 
with positive {(Tx) (blue region) and vice versa for neg- 
ative values (red region in Fig. [2];a) (iv)). This obser- 
vation can be understood, when rewriting the Bell state 
1^) = [(|0) + \l))\9x) + (|1) - |0))|e:c)]/2 in the eigenbasis 
{\gx)A^x)} of the measurement observable (Jx'- We note 
that for the field component (|0) + |1)) we find (X) > 
which is correlated with the state \gx)^ while le^^) is cor- 
related with (|0) - |1)) for which (X) < 0. Equivalently, 
the {(Jy) component is correlated with the sign of the P 
quadrature measurement as shown in Fig. [2|a) (v). 

Already in the raw measurement data we clearly ob- 
serve the expected qubit-photon correlations including 
their phase coherence. In order to further quantify the 



properties of the prepared Bell state we evaluate the sta- 
tistical moments {{a^)^a^(Ji) from the measured set of 
3-dimensional histograms using the methods presented in 
Refs. [^^. The resulting measured expectation values 
(colored bars) of products between the Pauli operators (Ji 
and photon field operators a, are compared with the 
theoretical values of an ideal Bell state (wireframes) up 
to order n+m = 4 in Fig.[2|b). Here we note that in com- 
parison to earlier measurements [I3j [15] , the increase in 
detection efficiency enabled by the parametric amplifier 
is essential for the measurement of higher order expecta- 
tion values which now also include products of qubit and 
photon field operators. 

The measured zeroth order moments {cFi) represent the 
Bloch vector of the qubit. Since all values are close to 
zero the qubit is, as expected, in the maximally mixed 
state when the photon part of |?/^) is traced out. The 
small finite value of (a^) is due to qubit decay during 
the time between state preparation and qubit tomogra- 
phy required for performing photon tomography in the 
same mode. Finite first order expectation values (a<j^) 
are due to the expected correlations between the equa- 
torial component of the Bloch vector and the phase of 
the photon field, already observed in the raw measure- 
ment data (Fig. [2|a) (iv)-(v)). The finite second order 
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FIG. 3: Density matrix for a two-photon entangled state. 
Real part of the measured (solid) and the ideal (wireframe) 
density matrix for the state |0) = |(|1) + |2))|^) + |(|1) - 
|2))|e) generated and measured with fidelity F — 80%. 

moments show that the single excitation is shared among 
qubit and photon field. Since the mean product of exci- 
tations {a^acjz) is close to the mean photon number {a^a) 
we find that whenever a photon is detected, the qubit is 
in the ground state for which Gz takes the value 1. We 
also find that all higher order moments with n-\-m = 3, 4 
are close to zero within their statistical errors, indicating 
that the photon field is a superposition of vacuum and 
single-photon states only [20 . In particular the mea- 
sured anti-bunching {{{a^)'^a^) = 0.023 ± 0.008) shows 
that there are no contributions of higher photon number 
states [13 . Moments of higher order can also be deter- 
mined from the measured histogram data (not shown), 
albeit with statistical errors which depend exponentially 
on increasing order [29 . 

We have also evaluated the density matrix p of the 
joint qubit-photon state from the measurement data us- 
ing a direct linear mapping from the moments to the 
density matrix elements [20 . In order to make use 
of the full measurement data and to guarantee a com- 
pletely positive density matrix we additionally apply a 
maximum-likelihood procedure which estimates the most 
likely density matrix from measured moments and their 
respective standard deviations up to order n + m = 8. 
This allows for reconstructing the density matrix in a 
10-dimensional Hilbert space including photon number 
states up to |n) = |4). As already expected from the 
vanishing fourth order moments, we find all number state 
populations with n > 1 close to zero (Fig.|2|c)). The co- 
herent superposition of the two contributing basis states 
|0e) and \lg) is reflected in the large off-diagonal ele- 
ments. Adjusting the overall local oscillator phase, the 
elements of the imaginary part of the density matrix (not 
shown) have been minimized to less than 0.023. The to- 
tal fidelity of the reconstructed state compared to the 
ideal Bell state is F = = 83%. The loss of fi- 

delity is dominantly due to qubit decay and decoherence 
during the 60 ns period between the state preparation 



pulse and the final tomography pulse, which determines 
the time at which the qubit state is characterized. As 
a measure of entanglement between the photon states 
{|0), |1)} and the qubit states {|^), |e)} we determine the 
concurrence C{p) = 0.72 which is clearly above the en- 
tanglement threshold C{p) = 0. 

To demonstrate the versatility of our state prepara- 
tion, detection and reconstruction scheme beyond exist- 
ing experiments we have prepared entangled states be- 
tween stationary qubits and multiple propagating pho- 
tons such as 10) = |(|1) + \2))\g) + ^(|1) - |2))|e). Mak- 
ing use of the third energy level / of the transmon [30] . 
we start the preparation sequence with a two-photon tt- 
pulse on the g-f transition which transforms the state 
from |0^) |0/). We then apply a sequence of two flux 
pulses to the transmon. Tuning the e-f transition in res- 
onance with cjr/27r during the first pulse, one excitation 
is swapped into the resonator |0/) |le). The second 
pulse entangles the state |le) with \2g) by tuning the g-e 
transition in resonance with the cavity for an appropriate 
time. An additional 7r/2 pulse applied at the g-e transi- 
tion frequency creates the state 10). Note that this state 
preparation sequence can be interpreted as the genera- 
tion of the separable state :^(|1) + |2)) (|^) + |e)) and 
an entangling controlled phase gate, which changes the 
sign of the two photon component |2) only if the qubit is 
in the excited state. The entanglement between photon 
field and qubit thus becomes apparent in the negative 
sign of the |2e) component in |0). We characterize the 
prepared state using the methods described above, which 
results in a final density matrix with fidelity F = 80% 
compared to the ideal one (Fig. [s]). The sign change 
apparent in the reconstructed state as the six negative 
(blue) elements clearly demonstrates the entanglement of 
the propagating multi-photon state with the stationary 
qubit. 

In our experiments, we have demonstrated the 
generation and detection of entanglement between a su- 
perconducting qubit and a propagating microwave field. 
The development of sensitive detection techniques for 
the measurement of photon-qubit quantum correlations 
is an important step towards using itinerant microwave 
photons as a quantum information carrier, e.g. for 
connecting spatially separated superconducting circuits 
or other systems interacting with microwave photons. 
The development of microwave photon counters ^ and 
the high level of control achievable over superconducting 
circuits put the realization of microwave photon based 
quantum network experiments within reach. 
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